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INTRODUCTION 


A  considerable  amount  of  the  energy  released  from  the  detonation  of  an  explosive  charge  can 
be  captured  using  a  metal  liner,  which  then  can  be  focused  to  form  a  projectile  for  boring  a  hole  in  a 
target.  A  typical  device,  referred  to  as  a  shaped  charge  (SC),  consists  of  an  explosive  charge 
formed  with  a  cavity,  which  is  lined  with  a  thin  layer  of  metal.  When  the  charge  is  detonated,  the 
resulting  pressure  collapses  the  metal  liner  and  forms  it  into  a  high-speed  projectile.  Depending  on 
the  explosive  and  the  mass  and  geometry  of  the  liner,  projectiles  can  attain  velocities  varying  be¬ 
tween  1.5  and  12  km/s.  When  such  a  high-speed  projectile  interacts  with  a  solid,  it  is  capable  of 
penetrating  it,  and  this  effect  is  employed  in  a  number  of  military  and  industrial  applications. 

Typically,  SCs  are  employed  either  (1)  at  short  (a  few  charge  diameters)  or  (2)  at  long  dis¬ 
tances  (tens  or  hundreds  of  charge  diameters)  from  the  target,  in  which  case  they  are  often  referred 
to  at  explosively  formed  penetrators  (EFPs).  Although  the  term  “explosively  formed  penetrators” 
implies  a  reference  to  a  projectile  (employed  as  a  penetrator)  produced  from  an  explosive  charge 
(i.e.,  shaped  charge),  an  accepted  usage  of  this  term  applies  to  a  charge  with  a  smooth  shallow  liner 
profile.  The  terminology  used  in  this  paper  will  follow  this  practice. 

Because  of  the  different  standoffs  for  different  applications,  times  available  for  the  projectile 
formation  can  be  markedly  different.  However,  in  order  to  maximize  target  penetration,  penetrator 
length  has  to  be  maximized.  Thus  for  most  applications,  the  given  standoff  (or  time)  to  form  a  pene¬ 
trator  before  striking  the  target  is  the  most  important  constraint  for  maximizing  projectile  length. 
These  considerations  determine  SC  (or  EFP)  effectiveness,  and  therefore,  play  a  crucial  role  in 
choosing  an  appropriate  liner  design  rationale. 

In  principle,  independent  of  the  liner  shape,  liner  acceleration  is  carried  out  very  rapidly  com¬ 
pared  with  the  time  allowed  for  projectile  formation;  therefore,  the  effectiveness  of  the  energy  trans¬ 
fer  from  the  explosive  to  liner  (or  resulting  projectile)  is  virtually  independent  of  the  target  standoff. 
Another  fact  that  is  important  in  liner  design  is  that  the  momentum  transferred  to  a  liner  is  a  fixed 
function  of  the  explosive  type,  liner  thickness  (mass),  and  the  charge  boundary  conditions,  and  the 
magnitude  of  the  momentum  is  not  significantly  affected  by  small  variations  in  liner  shape.  Never¬ 
theless,  the  dynamics  of  liner  collapse  and  projectile  formation  is  very  sensitive  to  the  liner  shape 
and  thickness.  Therefore,  the  adjustment  of  liner  shape  is  a  common  design  technique  for  the  con¬ 
trol  of  projectile  formation  time. 

For  example,  in  short-range  SC  applications,  where  the  time  for  the  formation  of  a  long  projec¬ 
tile  is  limited,  the  liner  is  designed  in  a  way  that  the  tip  of  the  resulting  projectile  (often  called  a  jet) 
moves  much  faster  than  the  tail.  This  results  in  a  stretching  projectile,  which  continues  to  gain 
length  even  as  its  tip  penetrates  the  target.  If  the  stretching  continues,  as  it  would  in  long  standoff 
applications,  the  projectile  would  ultimately  fracture.  This  can  significantly  degrade  target  penetra¬ 
tion.  Therefore,  for  long-range  (EFP)  applications,  the  design  criterion  is  to  achieve  long  and  stable 
(non-stretching)  projectile. 

The  employment  of  various  SC  designs  depends  on  specific  applications,  and  the  examples  of 
these  applications  are  numerous.  During  the  Second  World  War,  a  number  of  antitank  weapons 
(e.g.,  the  rocket-propelled  Bazooka)  were  developed  which  made  use  of  short-range  SCs  to  perfo¬ 
rate  thick  metal  plates  (ref.  1).  A  historical  outline  of  the  developments  in  shaped  charge  technology 
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is  given  by  Kennedy  (ref.  2).  Recent  long-range  EFP  developments  for  antitank  weapons  are  re¬ 
viewed  by  Homemann  et  al.  (ref.  3)  and  Weimann  (ref.  4).  Charges  lined  with  linear  wedges  (linear 
SC)  have  been  used  for  cutting  cables  and  bridge  beams  and  in  salvage  work  on  sunken  ships. 

Various  explosive  charges  are  commonly  used  in  mining  operations  to  break  large  rocks. 

Short  range  SCs  have  frequently  been  used  for  the  rapid  boring  of  holes  for  demolition  work.  After  a 
hole  has  been  produced  with  a  SC,  the  cavity  can  be  filled  with  explosive  for  further  blasting.  Some 
of  these  applications  (ref.  5)  have  stringent  requirements  for  the  bored  holes,  which  have  to  permit 
insertion  or  passage  of  a  fixed  diameter  main  demolition  charge.  This  paper  discusses  a  number  of 
the  issues  involved  in  the  development  of  a  special  single  SC,  which  produces  a  complex  penetrator 
that  should  be  ideally  suited  for  these  applications.  The  proposed  complex  penetrator  introduced  in 
this  work  has  two  separate  parts: 

•  A  stretching  central  solid  rod  penetrator  (referred  to  in  this  work  as  a  hole  “driller”) 

•  A  concentrate,  radially  stable  and  non-stretching  hollow  cylindrical  penetrator 
(referred  to  in  this  work  as  a  hole  “cleaner”) 

Our  experiments  confirmed  that  it  is  possible  to  produce  a  robust  penetrator  of  this  complex 
shape.  The  following  is  a  detailed  discussion  of  a  number  of  issues  involved  in  designing  this  SC. 


PENETRATOR’S  SHAPE 

Structures  that  are  encountered  in  demolition  operations  are  made  from  a  variety  of  geologic 
and  man-made  materials  and  are  usually  reinforced  with  steel.  When  a  simple  rod  penetrates  such 
a  target,  the  resulting  hole  diameters  vary  as  the  penetrator  penetrates  different  materials.  Ductile 
metals  are  ideally  suited  for  witnessing  this  effect.  Birkhoff  et  al.  (ref.  1),  using  a  stack  of  alternating 
steel  and  lead  plates,  demonstrated  a  striking  difference  in  hole  diameters  of  adjacent  plates:  hole 
diameters  in  the  lead,  which  is  extremely  weak,  were  four  to  five  time  larger  than  the  much  stronger 
steel.  Similarly,  when  the  target  material  is  inhomogeneous  and  consists  of  a  matrix  of  weak  brittle 
materials  containing  metal  inclusions,  these  effects  are  also  significant.  Appreciable  differences  in 
hole  diameters  between  the  embedded  steel  plates  and  the  concrete  were  reported  from  penetration 
experiments  with  gun  launched  penetrators  against  simulated  reinforced  concrete  targets  (ref.  6). 

In  this  type  of  target,  reinforcing  steel  occupies  only  a  small  fraction  of  the  total  hole  volume, 
but  nevertheless  controls  the  passage  of  the  demolition  charge.  This  requires  special  penetrator 
solutions  to  achieve  satisfactory  bore  holes.  When  a  hole  is  created  with  a  simple  rod  penetrator  in 
order  to  permit  passage  of  a  fixed  diameter  demolition  charge,  the  penetratoris  diameter  has  to  be 
increased.  This  severely  enlarges  the  hole  in  the  weaker  portions  of  the  target.  Such  holes  not  only 
witness  the  poor  efficiency  of  the  penetrator,  but,  most  importantly,  degrade  the  effectiveness  of  the 
demolition  charge  by  causing  excessive  venting  of  explosion  products.  Therefore,  to  ensure  proper 
performance  of  the  demolition  charge,  the  hole  has  to  be  created  with  two  penetrators:  (1)  a  hole 
driller,  which  creates  a  nominal  hole  in  the  weaker  material  and  (2)  a  hole  cleaner,  which  cuts  out 
remnants  from  the  harder  inclusions  (steel,  etc.).  A  variety  of  practical  approaches  are  possible,  and 
one  of  the  solutions  is  to  generate  these  penetrators  using  a  single  charge. 
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The  objective  here  was  to  investigate  a  single  charge  that  generates  two  penetrators:  a  central 
solid  rod  driller  and  a  concentric  hollow  cylinder  cleaner.  The  function  of  the  hole  driller  is  to  create 
a  nominal  diameter  hole.  The  deeper  the  burial  of  the  demolition  charge,  the  higher  is  the  degree  of 
the  confinement  of  the  explosive  and  the  resulting  transfer  of  the  momentum  to  the  target,  which 
ultimately  causes  it  failure.  Therefore,  increasing  hole  depth  increases  the  effectiveness  of  the 
demolition  charge.  Because  of  practical  considerations,  the  charge  standoff  is  limited  to  a  few 
charge  diameters.  Thus,  in  order  to  create  a  deep  hole,  the  hole  driller  has  to  stretch  as  it  pene¬ 
trates  the  target.  The  second  part  of  the  penetrator,  the  concentric  hollow  cylinder  acts  as  a  hole 
cleaner.  It  should  not  attempt  to  cut  through  the  main  target  material,  but  instead  should  follow  the 
hole  driller.  Therefore,  the  hole  cleaner  has  to  have  a  slower  velocity  than  the  hole  driller.  Since  a 
fixed  diameter  hole  is  required,  the  hole  cleaner  has  to  be  radially  stable;  i.e.,  it  should  not  continue 
to  collapse  or  expand  (radially)  as  it  moves  in  the  target. 

To  obtain  these  two  radically  different  penetrators  from  a  single  charge  is  not  a  simple  task. 
What  is  the  liner  shape  that  has  to  be  used  to  obtain  such  a  complex  projectile?  Figure  1  shows 
results  of  the  numerical  simulations  of  an  initially  attempted  design.  Both  parts  of  the  penetrator  are 
unsatisfactory.  In  fact,  the  simulations  suggest  that  the  peripheral  portion  not  only  fails  to  form  a 
robust  penetrator,  it  does  not  even  indicate  a  tendency  to  separate  from  the  central  part. 


NUMERICAL  SIMULATIONS 

To  produce  this  complex  penetrator,  the  principle  design  rationale  is  to  divide  a  liner  into  dis¬ 
tinct  central  and  peripheral  portions  that  are  to  generate  two  separate  penetrators;  hole  driller  and 
hole  cleaner,  respectively.  Typical  SC  and  EFP  liners  produce  penetrators  that  are  somewhat  simi¬ 
lar  to  the  two  required;  thus,  the  liner  profile  can  be  a  blend  of  SC  and  EFP  designs.  To  avoid  costly 
experimentation  with  different  liner  profiles,  numerical  simulations  are  instrumental  for  examination 
of  the  correlation  between  the  liner  shape  and  the  resulting  penetrator.  The  principle  step  in  this 
development  process  is  to  identify  basic  liner  parameters,  which  control  the  shape  and  robustness  of 
the  desired  penetrator.  Then,  based  on  observations  learned  from  the  simulations,  an  experimental 
charge  is  to  be  fabricated  and  tested,  and  the  experimental  data  are  to  be  compared  with  (and, 
hopefully,  prove)  the  tendencies  observed  in  the  analysis.  The  specific  objectives  of  numerical 
simulations  conducted  in  this  work  were  as  follows: 

♦  To  obtain  the  desired  hole  driller  with  the  required  penetration  parameter  such  as  tip 
velocity,  length,  stretching  rate,  diameter,  etc. 

♦  To  achieve  liner  separation  into  two  robust  parts  with  the  minimum  fragmentation 
losses 

♦  To  examine  liner  parameters  that  affect  formation  and  stability  of  the  desired  hole 
cleaner 

Numerical  Model 

The  dynamics  of  a  stretching  SC  liner  collapse  is  radically  different  from  a  stable  EFP  forma¬ 
tion.  The  projectile  formation  of  a  stable  EFP  is  accomplished  through  a  relatively  moderate  plastic 
deformation  of  the  liner,  similar  to  these  experienced  in  metal  forging.  Modeling  of  these  EFPs  is 
usually  carried  out  in  a  Lagrangian  reference  frame.  For  stretching  SCs,  the  plastic  deformation  of 
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the  liner  is  extremely  severe.  When  a  Lagrangian  reference  frame  is  employed  in  modeling  stretch¬ 
ing  SCs,  the  computational  mesh  is  severely  distorted.  This  has  an  adverse  effect  on  the  integration 
time  step  on  accuracy,  which  is  directly  related  to  the  grid  size  through  the  Courant  stability  condi¬ 
tion.  Therefore,  commonly,  SC  modeling  is  carried  out  using  Eulerian  computational  mesh,  which  is 
fixed  in  space  and  is  not  subject  to  the  Lagrangian  mesh  distortion  problems. 

In  our  problem,  the  penetrator  consists  of  both  a  stable  and  a  stretching  part.  Because  of  the 
mesh  distortion  of  the  stretching  portion,  purely  Lagrangian  modeling  would  not  be  acceptable. 
Although  using  a  Eulerian  reference  frame  would  solve  the  mesh  problems,  the  calculation  of  the 
lengthy  formation  of  the  stable  part  would  require  significant  computational  resources.  The  numeri¬ 
cal  simulations  presented  in  this  work  were  performed  using  the  CALE  (ref.  7)  computer  program. 

CALE  is  a  plane  two-dimensional  and  three-dimensional  axisymmetric  hydrodynamics  code 
basd  on  an  Arbitrary  Lagrangian  Eulerian  Formulation  of  the  governing  equations.  It  can  run  in  a 
pure  Lagrangian  mode,  a  pure  Eulerian  mode,  or  in  arbitrary  hybrid  mode.  CALE  is  capable  of 
modeling  multi-material  flows  and  allows  for  discontinuous  velocities  at  material  interfaces.  In  the 
Eulerian  reference  frame,  the  description  of  the  axisymmetric  kinematic  deformation  of  the  contin¬ 
uum  is  given  by  the  following  conservation  laws: 

Conservation  of  mass 


dp  dp  dp 
- h  Vr - H  Vz - H  p\ 

dt  r  dr  z  dz  1 


^  1  <9  /  \  dvz 

- (rvr)+ — - 

r  dr  d 
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z  J 


Conservation  of  linear  momentum 
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dr  *  dz  p\  dr  dz 
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Conservation  of  energy 


dE 


dE 
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dt  r  dr  z 


dE=  J_ 
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r  dr  dz  1 
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Body  forces  (such  as  gravity),  thermal  energy  sources,  and  heat  conduction  are  ignored. 

CALE’s  procedure  for  calculating  elastic-plastic  flow  is  based  on  the  methodology  developed 
by  Wilkins  (refs.  8  and  9).  The  stress  behavior  of  a  material  can  be  thought  of  being  composed  of  a 
stress  associated  with  uniform  hydrostatic  pressure  p ,  referred  to  as  the  spherical  part  of  a  stress 
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tensor,  plus  the  stress  associated  with  the  resistance  of  the  material  to  shear  distortion,  referred  to 
as  a  deviator  Sy .  In  describing  yielding  and  plastic  flow,  the  stress  contributions  that  are  due  to 

shear  distortion  are  limited  by  the  constitutive  equation  for  the  material.  The  stress  tensor  is  given 
then  by 


=~P+S,J 

The  constitutive  equation  for  the  liner  material  adopted  here  is  that  of  Steinberg-Guinan  (ref. 
10),  which  specifies  the  shear  modulus  and  the  yield  strength  as  a  function  of  pressure,  tempera¬ 
ture,  and  equivalent  plastic  strain.  The  hydrodynamic  response  of  liner  material  was  modeled  with 
the  standard  linear  polynomial  approximation,  usually  employed  for  metals.  The  liner  was  modeled 
with  a  standard  set  of  parameters  for  oxygen-free  high  conductivity  (OFHC)  copper,  which  were 
taken  from  Tipton  (ref.  9).  Details  for  these  models  and  their  implementation  can  be  found  in  Tipton 
(ref.  8). 

The  explosive  was  modeled  as  an  ideal  hydrodynamic  materia!  (s^  =  0),  serving  as  an  energy 

source.  The  procedure  for  modeling  the  explosive’s  energy  release  rate  used  the  Jones-Wilkins-Lee 
(JWL)  equation  of  state  and  a  burn  fraction  algorithm  based  on  a  constant  detonation  velocity  D0 
and  a  uniform  chemical  energy  density  e0 .  Parameter  values  for  LX-14  from  Tipton  (ref.  9)  were 
used  to  model  the  explosive. 

In  addition  to  the  specification  of  the  material  models,  problem  geometry  and  initial  and  bound¬ 
ary  conditions  have  to  be  specified  before  the  solution  procedure  can  be  initiated.  The  configuration 
of  an  axially  symmetrical  geometric  domain  in  r-z  space  is  shown  in  figure  2.  Domain  AXEXE3E2 
represents  the  explosive  booster  pellet,  E2E4E5E6E7E$ ,  is  the  main  explosive,  E^E7E6LXL2L3  is  the 
copper  liner,  RXR2R3R4 is  the  steel  confinement  ring,  and  the  rest  of  the  domain  AXA2A3A4  is  filled 
with  air. 

The  following  boundary  conditions  were  applied  along  the  boundary  of  domain  AXA2A3A4 


vz  =  0  along  A3A4 

cr„  =  0  along  AXA2  and  A2A3 

vr  =  0  along  AXA4  (line  of  axial  symmetry) 

CALE  is  capable  of  modeling  discontinuities  in  the  tangential  component  of  the  velocity  across  mate¬ 
rial  interfaces,  and  a  so-called  “slide  line  E6E7E%  between  the  explosive  and  the  liner  was  included 
in  the  model. 

The  explosive  was  detonated  at  the  point  Ax .  The  initial  size  of  the  computational  domain  was 
95  zones  in  the  z  direction  and  84  zones  in  the  r  direction.  With  the  explosion,  the  boundary  of  the 
spatial  domain  expands.  If  the  size  of  the  computational  domain  (i.e.,  number  of  zones)  was  to  be 
left  unchanged  with  time,  the  motion  of  the  stretching  projectile  that  reached  the  boundary  would 
have  to  be  resolved  on  a  limited  region  of  the  mesh.  This  would  introduce  great  inaccuracies  in  the 
calculations.  The  axial  component  of  the  velocity  at  the  boundary  A3A4  (the  air)  was  set  to  zero, 
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vz  =  0 ,  until  the  projectile  approached  it  at  a  distance  of  approximately  3  cm.  At  that  time,  3  cm  of 
new  mesh  containing  10  zones  (in  the  z  direction)  were  added  to  that  boundary.  Then,  as  the 
projectile  moved  onto  the  newly  formed  mesh  approaching  the  new  boundary,  the  procedure  was 
repeated.  The  periodic  addition  of  computational  mesh  was  performed  during  the  entire  course  of 
calculations. 

Hole  Driller 

When  the  explosive  detonates,  it  undergoes  vigorous  chemical  decomposition  that  generates 
pressures  in  the  order  of  two  to  four  hundreds  of  Kbars.  Since  the  ability  of  the  detonation  products 
to  transmit  shear  is  very  small,  the  liner  can  be  thought  to  be  accelerated  only  under  action  of  forces 
normal  to  its  surface.  By  changing  the  orientation  of  the  liner  surface  adjacent  to  the  explosive,  one 
can  change  the  relative  amounts  of  axial  and  radial  components  of  the  momentum  imparted  to  the 
liner.  Since  the  ensuing  motion  of  the  liner’s  material  depends  on  the  initial  amounts  of  axial  and 
radial  momenta,  the  adjustment  of  liner  angle  a  (fig.  3a)  is  a  common  technique  for  control  of  liner 
collapse  and  projectile  formation. 

Numerical  simulations  of  liner  collapse  for  different  liner  geometries  are  shown  in  figure  4.  In 
order  to  produce  the  required  stretching  hole  driller,  we  decreased  (as  compared  to  a  liner  shown  in 
fig.  1)  the  liner  angle  a .  Decreases  in  a  increase  the  radial  component  of  the  momentum  trans¬ 
ferred  to  the  liner.  This  causes  the  liner  to  undergo  more  intense  radial  deformation  (fig.  4).  Shaped 
charges  shown  in  figure  4  have  elliptic  liner  profiles,  and  their  liner  collapse  mechanisms  are  very 
similar  to  the  conical  SCs  (ref.  1).  In  order  to  illustrate  the  process  of  this  type  of  projectile  forma¬ 
tion,  one  can  apply  a  conservation  of  mass  principle  to  the  control  volume  occupied  by  the  liner. 
Liner  radial  contraction  generates  an  intense  flow  (or  jet)  in  the  axial  direction.  With  respect  to  the 
center  of  mass,  the  axial  flow  is  in  two  opposing  directions  and  the  difference  in  velocities  is  large. 
This  is  why  this  type  of  liner  collapse  mechanism  produces  a  stretching  solid  rod-like  projectile,  and 
when  the  projectile  (or  jet)  is  viewed  from  a  conventional  stationary  reference  frame,  axial  velocity  of 
the  tip  is  much  higher  than  of  the  tail. 

Liner  Separation  into  Two  Robust  Parts 

The  liner  angle  a  can  also  be  used  to  guide  the  direction  of  liner  collapse  in  such  a  way  that 
the  liner  would  separate  into  two  parts,  central  and  peripheral  (figs.  1  and  4).  This  is  accomplished 
with  the  following  technique.  Along  the  line  where  the  liner  is  to  break,  an  edge  is  formed.  At  the 
edge,  the  two  surfaces  are  given  significantly  different  collapse  angles  (if  the  center  has  to  stretch, 
a!  ox").  The  magnitudes  of  the  momenta  A/'andM'  imparted  to  the  two  parts  are  the  same,  but 
their  directions  are  different  (fig.  3b  and  c).  In  the  vicinity  of  the  edge,  when  the  two  momenta  are  in 
opposing  direction,  the  intense  plastic  flow  leads  to  “necking.”  This  ultimately  causes  the  liner  to 
separate  into  two  parts. 

What  liner  parameters  can  be  adjusted  to  control  the  separation  process?  Suppose  the  liner 
would  break  along  an  idealized  necking  plane  given  by  an  angle  p .  The  necking  is  driven  by  two 

forces  F  =  dM'/dt  and  F"  =  dM"/dt ,  and  the  separation  occurs  only  when  the  forces  pull  the  two 
parts  in  opposing  directions.  Therefore,  the  necessary  condition  for  the  separation  is  the 
(F1  •  n\F”  •«)<(),  which  requires  cc' </3<  a".  How  can  one  control  p?  By  changing  liner  thickness 
at  the  curvature  discontinuity.  As  a  first  approximation,  the  changes  in  the  thicknesses  do  not 
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change  the  forces  applied  to  the  liner.  However,  they  do  significantly  affect  the  relative  velocity  of 
the  two  liner  parts,  which  produces  the  plastic  flow  in  this  region.  Thus,  by  changing  the  relative 
thicknesses  of  the  two  parts,  one  can  control  the  orientation  of  the  necking  plane  ft . 

What  part  should  be  thicker?  In  order  to  increase  the  hollow  penetrator  robustness  the  ten¬ 
dency  should  be  to  increase  the  thickness  of  the  peripheral  part.  In  addition,  increases  in  thickness 
proportionally  decrease  penetrator  velocity,  which  is  also  an  important  factor  for  the  hole  cleaner 
application.  Thus,  the  liner's  peripheral  part  was  designed  to  be  thicker  than  the  central  part. 

The  employment  of  this  technique  alone  can  break  the  liner,  but  it  does  not  assure  that  the 
peripheral  part  would  produce  a  robust  penetrator.  The  following  considerations  were  employed  in 
the  design  of  the  robust  penetrator  shown  in  figure  4b.  Near  the  explosive  charge  boundary,  the 
detonation  products  are  free  to  escape  into  free  air,  which  causes  the  detonation  pressure  to  be 
relieved  very  quickly.  When  this  occurs,  the  amount  of  momentum  transferred  to  the  liner  is  de¬ 
graded,  primarily,  at  the  liner’s  periphery.  This  has  an  adverse  effect  on  the  formation  of  the  periph¬ 
eral  penetrator  resulting  in  excessive  stretching  that  leads  to  the  formation  of  the  peripheral  penetra¬ 
tor  resulting  in  excessive  stretching  that  leads  to  the  formation  of  a  conoidal  debris  cloud  instead  of 
the  desired  robust  penetrator  (fig.  4a).  The  amount  of  momentum  transferred  to  the  liner  periphery 
can  be  increased  with  a  steel  confinement  ring.  The  confinement  ring  prevents  the  rapid  escape  of 
the  detonation  products,  which  enhances  the  momentum  transfer  to  the  liner's  periphery  and  contrib¬ 
utes  to  formation  of  a  robust  penetrator.  Weickert  and  Gallagher  (ref.  11)  investigated  a  number  of 
confinement  ring  parameters  including  length,  thickness,  and  position;  our  choice  of  these  parame¬ 
ters  was  very  similar  to  their  findings.  When  a  confinement  ring  was  added  to  a  previous  design  and 
a  proper  choice  of  the  liner  angles  a'  and  a"  was  made,  numerical  simulations  indicated  formation  of 
the  robust  penetrator  shown  in  figure  4b. 

Hole  Cleaner  and  Idealized  Two-Mass  System 

The  hole  cleaner  is  formed  from  liner  peripheral  portion  (fig.  4b)  via  complex  plastic  deforma¬ 
tion,  which  is  simultaneous  with  translation  along  the  axis  of  symmetry.  In  this  process,  the  liner 
response  is  rather  complex:  liner,  which  is  initially  a  shallow  elliptical  shell,  is  subjected  to  more  then 
90-deg  inversion  into  a  conoidal  shell  that  undergoes  gradual  changes  in  shape.  To  idealize  this 
complex  motion,  we  represented  the  liner  as  a  two-mass  system  shown  in  figure  5.  In  this  two-mass 
system,  masses  m]  and  m2  represent  the  total  liner  mass  (the  peripheral  portion)  and  are  placed  at 
liner  edges  with  liner  angles  and  a2 ,  respectively.  Viscous  dampers  dj  and  d2  (fig.  5b)  oppose 
radial  motion  of  masses  mx  and  m2 ,  and  model  liner  resistance  to  radial  collapse.  Since  during 
plastic  deformation  liner  can  stretch,  the  masses  are  thought  to  be  connected  through  a  viscous 
damper  d2 ,  which  permits  the  distance  between  the  masses  to  change. 

The  total  momentum  transferred  from  the  explosive  to  liner  can  be  represented  with  the 
masses  ml  and  m2  given  initial  momenta  M\  and  M2 ,  respectively.  The  ensuing  motion  of  this 
relatively  simple  system  is  governed  by  a  limited  number  of  parameters  (initial  momenta  Mj  and 
M2 ;  and  damper  viscosities);  therefore,  from  the  analogy  with  the  motion  of  this  idealized  system, 

one  can  examine  a  number  of  important  correlations  that  can  guide  liner  design.  For  liner  design 
problems,  in  which  liner  design  changes  are  limited  only  to  small  changes  in  liner  shape,  the  primary 
parameter  determining  projectile  shape  is  the  direction  of  the  impulse  delivered  from  explosive  to 
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liner.  Damper  viscosities  are  parameters  that  model  liner  strength  and  relate  the  deformation  resis¬ 
tive  forces  to  displacements  produced  due  to  liner  plastic  flow.  Presence  of  the  dampers  dx,  d2,  and 
d3  in  a  model  is  crucial  for  analysis  of  a  number  of  the  hole  cleaner  parameters  including  the  radial 
stability  and  the  length.  However,  in  a  liner  shape  design  problem  where  the  concern  is  with  the 
relationships  between  liner  shape  changes  and  the  shape  of  the  resulting  penetrator,  all  the  damp¬ 
ers  can  be  omitted  and  the  model  can  be  substantially  simplified.  This  approximation  is  justified 
because  for  small  changes  in  liner  shape,  the  changes  in  the  dissipative  forces  modeled  by  dampers 
dh  d2,  and  d3  are  small  and  can  be  neglected.  Therefore,  for  such  a  model,  viscous  dampers  dx 
and  d2  are  omitted  and  a  rigid  link  is  substituted  for  damper  d3  (fig.  5c).  Their  shape  design  prob¬ 
lem  can  then  be  examined  by  drawing  analogy  between  changes  in  liner  geometry  and  variations  in 
momenta  M\  and  M2  of  an  equivalent  two  mass  system. 

Effect  of  Small  Liner  Shape  Changes 

Consider  an  axisymmetric  liner  surface  in  contact  with  the  explosive,  which  is  given  by  a  para¬ 
metric  curve 


*(<?)=**(£)*!  +<p2(Z)R2+V'  i(£)A  +  ¥2(^2 

where  0  <£<  1  is  a  dimensionless  coordinate  along  the  curve  R{%) -  r{%]er  +  z(d;)ez ,  Rx  and  R2  are 
vectors  at  the  end  points  (i.e.,  at  the  liner  edges),  Rx  and  R2  are  the  tangents  at  these  points,  and 
<P2^),  ¥\ (#).  and  ¥2^)  are  arbitrary  blending  functions  such  as 

^(0)-l,^(l)  =  0 

<Pl  (0)  =  0,  <32(0=1 

<p\  (0) = <Pi  (1) = 4>2  (°) = ^(0  =  0 
^(0)=1,^(1)=0 
^2(0)  =  0,^2(0)  =  0 
¥\  (0)  =  ¥\  0)  =  ¥ 2  (°)  =  ¥2  (0  =  0 


Now  consider  a  liner  design  problem,  in  which  liner  shape  is  changed  in  small  increments. 
Consider  that  these  small  changes  in  liner  profile  are  carried  out  through  small  changes  in  edge 
axial  positions  z,  and  z2  (fig.  6a)  and/or  liner  angles  tanaj  =  rx  I  zx  (fig.  6b)  and  tan  a2=r2lz2,  while 
the  radial  positions  rx  and  r2  of  the  edges  are  fixed.  Further,  suppose  that  blending  functions  are 
such  that  the  resulting  liner  profile  is  a  smooth,  small  curvature  function.  Clearly,  when  zx ,  z2 , 

and/or  ax  or  a2  are  changed,  both  R{£)  and  tan  a(<^)  =  r(^)/(i(^))_  1  change.  Since  projectile  forma¬ 
tion  is  extremely  sensitive  to  liner  profile  changes,  the  extent  of  variation  SR  arising  from  small 
variations  &2,  Sax,  and  Sa2  must  be  discussed  thoroughly. 

First,  consider  the  magnitude  of  the  total  momentum  transferred  to  the  liner,  which  can  be 
expressed  as 


1  r 

=  2zrJ  J  p{%,  iy(£)d%dt 


00 
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where  t  is  the  time,  p(g,  t),  t  <  r;  p  =  0,  t  >  t  is  the  pressure  exerted  on  the  liner  by  the  explosive. 
Since  the  momentum  is  directly  proportional  to  the  radial  component  r(g)  of  liner  profile,  the  primary 
effects  of  liner  shape  changes  are  the  variations  Sr . 

Consider  first  that  the  effect  on  changes  and  Sz2  in  the  axial  positions  of  the  liner  edges 
(fig.  6a).  Since 


Ks)=  Pi  (&i  +  <Pl(&2  +  Y\  (&]  +  Yl(£)h 

clearly,  Sr(g,  5z\ ,  Sz2 )  =  0 .  Now  consider  a  situation  where  liner  profile  is  adjusted  by  changing  the 
liner  angle  (fig.  6b).  Since  a  small  change  5ocx  changes  only  the  direction  of  the  tangent  vector 
=||/?i||efi ,  the  resulting  variation  in  position  of  point  g  due  to  a  change  Sa\  is 


SR(g,Sa\)=n(Z)\R\\s(etX) 

Noting  that  the  rotational  change  of  the  unit  vector  en  is  given  by 

S(en)  =  Sa \eqxen 


the  variation  in  liner  profile  is  then 

SR{^, Sal)=y/l (g) || Saxeq  x en 

Figure  6c  illustrates  typical  shapes  of  the  blending  function  ^(£).  For  smooth,  small  curvature 
profiles  •  From  figure  6d,  it  can  be  shown  that 

Sr(^.  Saj )~  01(g)  l^i  ||  Sa'j  cos  ocq 

where  angle  a0  is  defined  in  figure  6d  and  arctan((r2  -  )l(z2  ~z\))<ao  <cc\-  For  typical  shallow  liner 
profiles,  cosa0  is  small.  Thus,  for  and  small  SahSr(g,Sai )  is  negligible.  Similarly,  it  can  be 

shown  that  Sr(g,Sa2)  is  also  small,  and  we  conclude  that  any  small  liner  shape  adjustments  &y , 

&2 .  .  and  Sa2  do  not  affect  r(g). 

Now  consider  the  changes  in  pressure  profile  p(g,t)  due  to  small  liner  shape  changes.  Since 
the  radial  variations  Sr  are  negligible,  a  change  in  pressure  at  the  interface  between  the  explosive 
and  the  liner  can  be  only  due  to  the  axial  variations  Sz .  A  typical  pressure  history  at  a  point  at  the 
interface  between  the  explosive  and  the  liner  can  be  generically  represented  by  a  curve  rapidly 
decreasing  off  a  Chapman-Jouguet  state.  The  profile  of  this  curve  reflects  the  pattern  of  interaction 
of  rarefaction  waves  originated  along  the  entire  charge  boundary,  including  the  liner  free  surface. 

With  a  change  in  the  geometry  of  a  boundary,  the  distances  traveled  by  the  rarefaction  waves 
change;  therefore,  any  liner  profile  modifications  are  expected  to  affect  the  profile  of  the  pressure 
history. 
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When  the  charge  diameter  is  fixed,  the  most  significant  change  in  the  boundary  is  in  the  axial 
position  of  the  liner  free  surface.  Undoubtedly,  this  changes  the  time  of  arrival  of  the  detonation 
wave  at  the  explosive  liner  interface.  For  example,  for  a  plane  detonation  wave  propagating  with  a 
velocity  £>0 ,  a  change  in  the  detonation  wave  time  of  arrival  is  given  by 


For  small  dz  and  a  shallow  liner  profile,  relative  changes  in  distances  between  pairs  of  arbitrary 
points  at  the  interface  between  the  explosive  and  the  liner  are  negligible.  Therefore,  it  is  reasonable 
to  assume  that  the  pressure  history  profile  is  not  significantly  affected  by  a  small  change  in  the  liner 
free  surface.  Hence,  the  pressure  for  a  modified  liner  can  be  expressed  as 

p{4,t+8r)*  p{$j) 

Similar  assumptions  can  be  applied  to  the  relaxation  time  t ,  which  is  also  expected  not  to  vary 
significantly  with  slight  changes  in  the  liner  geometry.  Thus,  the  value  of  the  integral 


\p{£,t)dt 

o 

will  be  unchanged.  Since  r{£)  is  not  affected  by  small  changes  in  liner  geometry,  we  can  conclude 
that  the  magnitude  of  the  total  momentum  M0  transferred  to  the  liner  is  invariant  to  slight  liner  modi¬ 
fications. 

The  radial  and  the  axial  components  of  the  total  momentum  are  given  as 

1  T 

Mqv  =  2/r J  J  p  cos  a{^)  r{£)d%dt 
00 

1  r 

M0z=  2xj  f  p  sin  a(£)  rfyd&t 
00 

With  adjustments  in  liner  shape,  the  liner  angle  a{£)  can  change  substantially,  which  effects  the 
ratio  between  the  axial  and  the  radial  components  of  the  momentum  transferred  to  the  liner.  Thus, 
the  primary  effect  of  small  liner  shape  adjustments  is  to  change  the  direction  of  the  total  momentum, 
while  its  magnitude  remains  the  same.  Similarly,  we  suggest  that  liner  shape  adjustment  carried  out 
by  changing  liner  angles  aj  and  a2  is  analogous  to  changing  the  “ratio”  between  the  momenta  Mx 
and  M2  of  the  two-mass  system,  while  their  total  magnitude  is  fixed.  Results  of  numerical  simula¬ 
tions  shown  in  figure  7  and  discussed  next  confirm  the  validity  of  this  analogy. 
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Numerical  Simulations  of  Effects  of  Small  Liner  Shape  Changes  on  Hole  Cleaner 

Refer  to  figure  7a  through  e,  which  shows  the  effects  of  small  liner  geometry  variations  on  hole 
cleaner  design,  all  calculated  at  130  ps.  All  liner  shape  adjustments  presented  here  were  achieved 
by  changing  peripheral  liner  thickness  and  liner  angles  a}  and  a2 ,  while  the  radial  positions  of  the 
outer  edge  and  the  central  portion  (that  produces  the  hole  driller)  were  fixed.  It  is  interesting  to  note 
that  all  calculations  give  essentially  identical  hole  driller  shapes.  This  illustrates  the  usefulness  of 
the  employed  charge  design  rationale  that  is  based  on  a  concept  of  mutual  independence  of  the 
formation  processes  of  the  hole  driller  and  the  hole  cleaner. 

Figure  7a  presents  calculations  wherein  liner  angle  a2  was  changed  by  5  deg,  which  produced 
approximately  16%  relative  radial  “translation”  of  the  hole  cleaner.  This  liner  profile  was  designed 
using  a  small  axial  shift  of  the  liners  outer  edge,  while  ax  and  liner  thickness  were  fixed.  This  type 
of  liner  shape  adjustment  can  be  expressed  as  a{g)±y ,  where  depending  on  the  shift  direction,  the 
liner  angle  is  either  increased  or  decreased  by  an  angle  y .  For  typical  liner  profiles  with 
60°  <«(£)<  90°  and  small  y 

sin(a±/)«sin(a) 

cos(a±/)«cos(a)+/ ,  60° <«(£)< 90°,  small  y 

with  a  maximum  of  12%  error  bound.  The  magnitudes  M0r  and  M0z  are  directly  proportional  to  the 
values  of  sm(a±y)  and  cos(«±/),  respectively,  therefore,  we  can  draw  a  conclusion  that  such 
adjustments  can  only  change  the  radial  component  of  the  momentum,  while  the  axial  component 
remains  approximately  the  same.  Thus,  the  radial  translation  effect  observed  in  figure  7a  is  due  to 
changes  in  the  radial  component  of  the  momentum,  which  caused  significant  change  in  the  hole 
cleaner  diameter.  In  terms  of  the  idealized  two-mass  system,  such  liner  shape  adjustment  is  equiva¬ 
lent  simply  to  a  proportional  change  of  both  MXr  and  M2r  while  ||A/1|+||M2||  is  fixed. 

From  penetrator  design  considerations,  it  can  be  desirable  to  decrease  the  hole  cleaner  veloc¬ 
ity  without  affecting  its  shape.  From  an  analogy  with  the  idealized  two-mass  system,  a  simple  solu¬ 
tion  is  to  proportionally  change  both  masses  mx  and  m2 ;  in  a  liner  design  this  can  be  readily  imple¬ 
mented  by  changing  liner  thickness.  Figure  7b  compares  two  designs  with  identical  liner  profiles 
and  different  thicknesses.  Results  in  figure  7b  show  a  distinct  difference  in  the  axial  position  of  the 
two  hole  cleaners  produced  from  these  designs,  while  their  shape  is  essentially  the  same. 

Consider  another  hole  cleaner  design  problem  in  which  the  desired  shaped  is  a  hollow  cylin¬ 
der.  Using  the  two-mass  system  analogy,  a  poor  penetrator  shape  can  be  explained  from  the  mo¬ 
tion  of  the  idealized  masses  mx  and  m2 .  Motion  of  these  masses  can  be  viewed  as  a  “clockwise” 
rotation  about  their  center  of  mass  combined  with  the  center  of  mass  translation.  As  seen  in  figure 
7,  at  130  ps  the  penetrator  had  formed  in  a  way  that  its  wall  is  at  some  angle  with  the  z  axis.  In 
order  to  increase  this  angle  to  the  desired  180  deg,  the  liner  wall  rotation  rate  should  be  changed. 
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Let  us  now  consider  the  factors  affecting  the  liner  wall  rotation  rate.  First,  consider  the  ideal¬ 
ized  two-mass  system.  Conservation  of  angular  momentum  requires  that 

dH 

8  =  ~dT 

The  torque  G  that  is  initially  imparted  to  the  masses  can  be  expressed  as 

Go  -(1-0*2) 

where  F,  =  dMx  Idt,  F2  =  dM2 I  dt,  4cm  is  the  position  of  the  center  mass  of  the  system  of  masses  ml 
and  m2 ,  and  /  is  the  distance  between  the  masses,  whose  relative  orientation  is  given  by  a  unit 
vector  <?i .  After  the  initial  angular  impulse  is  delivered  to  the  system,  and  in  the  absence  of  dissipa¬ 
tive  forces,  G  =  0 .  Hence,  integrating  the  angular  momentum  balance  equation  yields  the  rotational 
rate 


lio  =  H(o) 

where  1  is  rotational  inertia  of  the  system.  The  angular  momentum  H(0)  at  /  =  0  can  be  expressed 
as 


(  * 

H(0)  =  le1x 

1 

1 

KT 

tj 

l  0 

0  ) 

or  simply 

H(0)  =  lei  x  i£cnM\  -  (l  -  4cm  )^2  ) 

Thus,  for  a  given  system  with  fixed  £cw ,  the  magnitude  of  the  angular  momentum  depends  on  the 
vector  difference  between  the  momenta  Mx  and  M2 .  Therefore,  the  system  rotation  rate  can  be 
varied  by  changing  relative  directions  of  momenta  M\  and  M2 .  Angular  momentum  transferred  to  a 
liner  is 

1  X 

H=2rc\\prg){Rg)-Rcm)xngWt 

00 

where  n(%)  is  a  unit  normal  of  liner  surface,  Rcm  is  the  position  of  the  liner  center  of  mass,  and 
R(^)=  r{£)er  +z(^)ez .  Let  us  assume  that  for  these  liner  design  problems,  small  changes  in  liner 
shape  R{%)  produce  negligible  changes  in  angular  inertia.  Thus,  similarly  to  the  two-mass  system, 
where  a  change  in  rotation  rate  is  achieved  by  changing  the  momenta  directions,  the  rate  of  the  liner 
wall  rotation  can  be  influenced  primarily  by  changing  the  liner  surface  orientation  nig). 
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In  design  E  shown  in  figure  7c,  the  shape  of  the  original  liner  C  was  changed  using  an  ap¬ 
proximately  5-deg  increment  («!)  and  decrement  (< i2 )  of  the  liner  edge  angles,  while  the  positions 

of  both  edges  were  fixed.  Note,  that  since  the  “average”  liner  angle  J/ (a\  +  a2)  remained  the  same, 

the  total  momentum  (linear  momentum)  transferred  to  the  liner  was  unchanged.  However,  due  to 
the  higher  liner  curvature  of  design  E ,  the  total  angular  momentum  transferred  to  the  liner  was 
increased.  Calculated  penetrator  shapes  produced  from  the  designs  E  and  D  are  compared  in 
figure  7c.  In  agreement  with  the  previous  discussion,  an  increase  of  the  angular  momentum  in¬ 
creased  the  liner  wall  rotation,  and  figure  7c  shows  that  the  wall  of  design  E  “straightened  up.”  In 
this  particular  design,  liner  D  is  lighter  than  E ,  therefore,  the  position  of  liner  D  is  slightly  shifted  to 
the  right.  Figure  7d  compares  the  design  E  with  design  F ,  which  has  the  same  thickness,  but 
substantially  different  curvatures  (due  to  an  approximately  1 1  deg  difference  between  the  corre¬ 
sponding  liner  edge  angles).  In  this  case,  the  difference  between  liner  curvatures  is  much  higher, 
and  the  effect  of  rotation  rates  is  even  more  pronounced. 

In  the  preceding  discussion,  we  omitted  all  viscous  dampers  from  the  two-mass  system.  Let 
us  first  consider  the  effect  of  the  viscous  damper  d2 ,  which  models  liner  stretching  and  permits  the 
distance  between  the  masses  to  change.  If  in  the  idealized  two-mass  system  the  rigid  link  between 
masses  is  released  (i.e.,  replaced  with  a  damper  with  zero  viscosity),  the  stretching  rate  can  be 
expressed  as 


Mj  •  <2]  M2  ■ 
m{  m2 


Since  the  rate  of  stretching  is  proportional  to  the  relative  difference  between  the  velocities  of  the 
masses,  changing  the  mass  ratio  effects  stretching.  In  the  designs  E  and  D ,  we  changed  liner 
thickness  uniformly  by  5%,  which  resulted  only  in  slight  changes  of  the  center  of  mass  velocity  (fig. 
7b).  Now,  consider  design  G  whose  liner  angle  profile  is  identical  to  E ,  but  the  thicknesses  are 
different:  at  edge  1  by  57%  and  at  edge  2  by  the  same  5%  (as  in  designs  E  and  D).  Calculations 
of  hole  cleaner  shapes  produced  from  these  designs  are  compared  in  figure  7e.  Calculations  show 
not  only  increases  in  the  center  of  mass  velocity,  but  in  addition,  a  significant  increase  in  length. 
Since,  in  general,  penetrator  effectiveness  is  enhanced  with  increased  in  length,  this  might  be  a 
useful  factor  to  consider  in  designing  hole  cleaners. 

A  requirement  that  is  important  for  the  effectiveness  of  a  hole  cleaner  is  its  radial  stability;  i.e., 
after  the  desired  hollow  cylinder  with  equal  front  and  rear  diameters  D\  and  D2 ,  respectively,  if 
formed,  radial  deformations  should  cease.  As  discussed,  formation  of  the  hole  cleaner  is  accom¬ 
plished  through  rotation  of  the  liner  wall,  and  liner  rotation  with  continue  until  the  angular  momentum 
initially  imparted  to  the  liner  is  dissipated  in  plastic  flow.  Several  suggestions  to  obtain  a  radially 
stable  design  can  be  offered.  Consider  that  the  design  objective  has  produced  a  hole  cleaner  with 
given  diameter  D  at  a  given  standoff  (or  formation  time).  The  first  step  will  be  to  design  a  liner  that 
produces  a  penetrator  with  front  diameter  Dj  «  D  at  which  time  the  liner  rotation  has  ceased.  In 
terms  of  the  two-mass  system,  such  a  design  can  be  obtained  by  proper  balance  between  the 
damper  viscosities  d{  and  d2 ,  and  the  initial  amount  of  the  angular  impulse  delivered  to  the  system. 
Since  the  damper  viscosities  represent  the  forces  resisting  plastic  deformation,  the  corresponding 
liner  parameters  are  the  liner  material  strength  and  thickness.  Suppose  that  such  a  design  is 
achieved,  but  D2>  E\»D.  Once  such  a  design  id  established,  the  second  step  is  to  try  increasing 
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the  hole  cleaner  wall  rotation  by  changing  liner  curvature.  This  is  quite  feasible,  since  in  the  result 
shown  in  figure  7d  more  than  a  45-deg  rotation  change  (with  the  rotation  still  continuing)  was 
achieved  by  an  1 1-deg  change  in  the  liner  edge  angles.  Changing  liner  curvature  might  also  disturb 
the  achieved  balance  between  angular  impulse  and  dissipation,  which  might  result  in  [\>  D .  When 
Z)j  >  D  occurs,  we  continue  iteration  from  the  first  design  step  and  try  to  decrease  E\  by  increasing 
liner  thickness.  To  achieve  a  suitable  design  will  definitely  require  more  than  one  design  iteration 
and  the  performance  of  each  iteration  should  be  verified  experimentally. 


EXPERIMENTS 

The  explosive  charge  was  cast  from  molten  Octol  70/30  explosive  using  a  sedimentation  tech¬ 
nique.  This  technique  allows  the  gravity  settling  of  heavier  HMX  particles  in  the  melted  explosive, 
which  permits  a  higher  concentration  of  more  energetic  HMX  particles  in  the  main  portion  of  the 
charge  (close  to  the  liner),  while  the  less  energetic  TNT  flows  upward  toward  the  riser. 

The  charge  casting  procedure  was  as  follows: 

The  device  was  assembled  from  the  aluminum  casing,  the  liner,  and  the  steel  confine¬ 
ment  ring,  which  were  bonded  together  using  an  adhesive.  The  Octol  70/30  explosive  was  melted 
and  heated  to  approximately  96°  to  104°C.  A  vacuum  about  686  mm  of  Mercury  was  applied  to  the 
melt  kettle  for  15  min.  The  explosive  was  poured  into  the  preheated  charge  assembly  and  vibrated 
for  10  min.  The  assembly  was  maintained  at  a  minimum  temperature  of  96°C  for  at  least  1.5  hrs. 
Compressed  air  at  ambient  temperature  was  then  applied  to  the  liner  surface  while  the  riser  was 
maintained  above  93°C.  This  allowed  the  explosive  to  solidify  from  the  liner  upward  toward  the  riser. 
The  heat  source  to  the  riser  was  later  removed,  and  the  charge  was  allowed  to  cool  slowly  to  ambi¬ 
ent  temperature.  The  loading  funnel  and  the  riser  were  then  removed  and  the  surface  of  the  explo¬ 
sive  was  hand  smoothed  to  final  design  dimensions.  The  1.9-cm  long  and  2.54-cm  diameter  CH-6 
booster  pellet  was  placed  into  a  cylindrical  cavity  (0.3-cm  deep)  machined  in  the  cast  explosive.  A 
Reynolds  RP-2  detonator  was  centered  on  the  booster  with  a  machined  detonator  holder. 

The  objective  of  this  experiment  was  to  prove  the  feasibility  of  producing  two  robust  penetra- 
tors  from  a  single  charge.  Therefore,  clear  images  of  the  penetrator’s  shape  during  projectile  forma¬ 
tion  (starting  from  a  relatively  early  time)  were  crucial  to  the  success  of  the  experiment.  The  follow¬ 
ing  considerations  motivated  the  usage  of  flash  (i.e.,  short  exposure)  radiography.  The  explosion 
generates  an  opaque  cloud  of  detonation  products  and  miscellaneous  debris,  which  blocks  the 
projectile’s  view  at  early  projectile  formation  times.  This  phenomenon  limits  the  utility  of  optical 
photography.  The  density  of  the  debris  cloud  is  low,  so  that  the  cloud  is  transparent  even  to  low 
intensity  x-ray  radiation.  This  makes  flash  radiography  ideally  suited  for  this  experiment.  As.  the 
penetrator  formed,  three  x-ray  tubes  exposed  x-ray  film,  which  recorded  the  projectile  images. 

A  schematic  of  the  experimental  setup  is  shown  in  figure  8.  The  equipment  was  conducted 
inside  an  explosion  proof  chamber.  The  three  x-ray  tubes  were  installed  in  special  ports  in  the  wall 
and  in  the  ceiling  of  the  chamber.  X-ray  tube  A  (fig.  8)  was  installed  in  the  ceiling,  directed  to  the  x- 
ray  film  cassette  A,  and  provided  the  top  view  of  the  penetrator.  The  x-ray  tubes  B  and  C  were 
installed  in  the  same  wall  port  and  provided  two  side  views  of  the  penetrator.  Tubes  B  and  C  were 
positioned  in  a  plane  orthogonal  to  the  projectile  path  and  at  an  angle  to  each  other  (tube  C  was  in  a 
horizontal  plane).  Tubes  B  and  C  were  pointed  to  the  same  cassette  B;  the  lower  portion  of  the 
cassette  recorded  the  image  from  tube  B;  while  the  upper  portion  -  from  tube  C. 
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The  x-ray  port  positions  are  fixed  in  the  explosion  chamber,  and  they  had  to  be  used  in  a  way 
that  provided  the  fullest  views  for  certain  x-ray  exposure  times.  On  the  other  hand,  the  desired 
exposure  times  have  to  represent  characteristic  stages  of  projectile  formation.  Thus,  x-ray  exposure 
times  of  105  ps  and  200  jj.s  were  selected  as  a  compromise  between  the  fixed  x-ray  port  positions 
and  the  desired  times  based  on  the  numerical  simulations.  These  times  were  used  to  trigger  x-ray 
tubes  A  and  B,  respectively.  The  third  x-ray  exposure  at  220  ^s  (for  tube  C)  was  set  a  short  time 
after  the  second  exposure,  which  was  necessary  for  the  reliable  measurement  of  projectile  veloci¬ 
ties. 


The  longitudinal  position  of  the  charge  in  the  chamber  was  selected  in  reference  to  the  first  x- 
ray  port  position  (tube  A)  and  was  based  on  the  projectile  velocities  calculated  in  the  numerical 
simulations.  The  charge  was  positioned  in  a  light  Styrofoam  stand.  A  stack  of  fifteen  25.4  x  25.4-cm 
square  5-cm  thick  steel  witness  plates  was  positioned  109  cm  from  the  edge  of  the  charge  on  an 
auxiliary  wood  frame,  which  also  supported  the  charge  stand.  A  7.5-cm  diameter  steel  calibrating 
rod,  with  fiducial  marks  5  cm  apart,  was  positioned  along  the  expected  projectile’s  path.  The  film 
was  then  irradiated  with  a  small  doze,  which  reproduced  the  pattern  of  fiducial  marks  on  the  film. 

This  produced  a  convenient  spatial  reference  frame,  which  simplified  data  reduction  procedures 
such  as  the  reference  of  the  image  positions,  the  x-ray  magnification  factors,  etc.  After  the  expo¬ 
sure,  the  calibrating  rod  was  removed.  During  the  experiment,  the  films  were  exposed  again  and  the 
projectile’s  images  were  superimposed  with  the  images  of  the  fiducial  marks. 

The  radiographic  images  of  the  penetrator  exposed  at  105  ps,  200  |xs,  and  220  p,s  are  shown 
in  figure  9.  All  radiographs  show  that  the  liner  separated  into  two  robust  parts.  At  105  jos,  the  site  of 
separation  is  still  well  preserved  on  the  central  portion  of  the  projectile.  The  characteristic  shape  of 
this  site  resulted  from  intensive  plastic  flow  (i.e.,  necking)  in  the  vicinity  of  the  liner  curvature  discon¬ 
tinuity,  and  it  agrees  very  well  with  the  numerical  simulations  (figs.  4  and  7).  The  later  radiographs 
of  200  (xs  and  220  jxs  indicate  formation  of  a  ring  of  debris  at  the  site;  i.e.,  the  material  which  formed 
that  neck  had  fractured. 

In  order  to  trace  the  evolution  of  the  projectile  shape  with  time,  four  characteristic  points  were 
identified  (fig.  9):  points  1  and  2  represent  an  idealized  shape  of  the  home  cleaner;  points  3  and  4 
are  the  respective  positions  of  the  tail  and  the  tip  of  the  hole  driller.  Tabulated  r-z *  coordinates  of 
these  points  are  given  in  figure  9.  The  evolution  of  penetrator  shape  agrees  very  well  with  the  re¬ 
sults  predicted  from  the  numerical  simulations.  Based  on  radiographs  at  200  ^s  and  220  p.s,  the  tip 
and  the  tail  of  hole  driller  mover  with  velocities  4.3  km/s  and  0.9  km/s,  respectively.  These  velocities 
agree  very  well  with  the  numerical  simulations. 

Based  on  the  data  in  figure  9,  the  average  axial  velocity  of  point  1  (the  front)  of  the  hole 
cleaner  between  105  jxs  and  200  |xs  was  2.93  km/s,  while  between  200  ps  and  220  ps  it  was  2.90 
km/s.  In  the  same  time,  the  axial  velocity  of  point  2  (the  rear  of  the  hole  cleaner)  increased  from 
2.55  km/s  to  2.70  km/s.  These  data  agree  very  well  with  the  tendencies  observed  in  the  numerical 
simulations:  formation  of  the  hole  cleaner  is  accomplished  through  a  deceleration  of  the  front  and 
acceleration  of  the  rear.  As  the  hole  cleaner  continues  to  form,  the  radial  velocities  of  its  rear  and 
front  change  from  -0.18  km/s  to  -0.10  km/s  and  from  -0.24  km/s  to  +0.05  km/s,  respectively.  The 
data  indicate  significant  rotation  of  the  hole  cleaner,  which  confirms  the  formation  mechanism  ob¬ 
served  in  numerical  simulations.  When  the  rotation  reaches  a  pint  when  the  diameter  of  the  pene¬ 
trator  front  starts  to  increase  (i.e.,  the  front  starts  to  stretch  circumferentially),  it  adversely  affects  the 
robustness  of  the  hole  cleaner.  Ultimately,  this  rotation  will  lead  to  the  fracture  of  the  hole  cleaner. 
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A  photograph  of  the  first  four  witness  plates  and  a  schematic  of  the  hole  profile  are  given  in 
figure  10,  which  shows  two  distinct  hole  portions  produced  by  the  hole  driller  and  by  the  hole 
cleaner.  Since  the  hole  produced  by  the  cleaner  is  well  rounded  and  symmetrical,  it  is  assumed  that 
the  penetrator  did  npt  fracture.  On  the  other  hand,  because  of  the  continuing  rotation,  the  perform¬ 
ance  of  the  hole  cleaner  was  degraded.  The  supporting  considerations  are  as  follows:  The  target 
was  located  approximately  40  cm  from  the  position  of  the  last  radiographic  image  (220  ps)  of  the 
hole  cleaner  and  at  this  time  its  length  was  6.7  cm.  Based  on  ideal  hydrodynamic  theory  of  penetra¬ 
tion  (ref.  1),  the  penetration  depth  can  be  approximately  estimated  as  P  =  L'lppl  pt,  where  L  is  the 

length  of  the  penetrator  and  pp  and  pt  are  the  densities  of  the  penetrator  and  the  target.  If  the 

hole  cleaner  was  stable,  simple  calculations  give  an  estimate  for  penetration  to  be  at  least  7.2  cm, 
which  is  approximately  40%  deeper  than  the  depth  witnessed  in  the  experiment.  Thus,  we  con¬ 
clude,  that  at  the  time  of  the  encounter  with  the  target  (»360  ps)  the  hole  cleaner  length  was  de¬ 
graded  to  approximately  4.7  cm.  The  reason  for  this  decrease  in  length  was  continuing  rotation  of 
the  hole  cleaner. 


CONCLUSIONS 

Rapid  hole  boring  for  demolition  operations  requires  a  complex  penetrator  consisting  of  two 
separate  parts:  a  hole  driller  and  hole  cleaner.  To  produce  this  complex  penetrator  from  a  single 
charge  is  not  a  simple  task.  We  provide  the  validity  of  a  liner  design  rationale  wherein  formation  of 
the  hole  driller  and  hole  cleaner  are  mutually  independent,  and  their  shape  can  be  adjusted  using 
incremental  changes  in  liner  profile.  We  identified  a  number  of  parameters  tat  are  crucial  for  pre¬ 
venting  liner  fracturing  during  separation  into  two  parts  and  ensuring  the  robustness  of  the  resulting 
penetrator.  The  experiments  proved  the  validity  of  the  approach  and  indicated  the  feasibility  of 
producing  a  robust  penetrator  of  this  complex  shape.  The  hole  driller  produced  in  the  experiments  is 
quite  satisfactory  for  the  hole  boring,  while  the  hole  cleaner  requires  additional  adjustment.  The 
deficiencies  of  the  achieved  hole  cleaner  design  (such  as  poor  shape  and  radial  instability)  are 
typical  for  many  interim  hole  cleaner  design  that  were  observed  in  numerical  simulations  and  dis¬ 
cussed  in  this  work.  Formation  of  the  hole  cleaner  is  accomplished  through  rotation  of  the  liner  wall, 
and  this  rotation  continues  until  the  angular  momentum  initially  imparted  to  liner  is  dissipated  in 
plastic  flow.  Since  the  initial  liner  angular  impulse  and  its  dissipation  can  be  controlled  by  adjust¬ 
ments  in  liner  profile  and  thickness,  a  suitable  design  may  be  achieved  via  proper  incremental  liner 
shape  changes.  This  may  be  accomplished  using  a  number  of  techniques  developed  in  this  work 
and  would  require  a  few  analytical  and  experimental  design  iterations. 
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Figure  1 

Numerical  simulations  of  an  initially  attempted  design 
Both  the  central  and  the  peripheral  parts  of  the  penetrator  are  unsatisfactory. 


Confinement  Ring 


Figure  2 

Initial  configuration  of  explosive  charge  and  liner 
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Figure  3 

The  adjustment  of  the  liner  angle  a  is  a  common  technique  for  control  of  liner 

collapse  and  projectile  formation 

Changing  liner  angles  a'  and  a ”  and  the  orientation  of  the  necking  plane  p ,  one  can  control  liner 
separation  into  two  pats.  The  liner  in  (b)  does  not  satisfy  the  condition  for  liner  separation 

(a'  <p<a"),  while  the  liner  in  (c)  does. 


Figure  4 

Numerical  simulations  of  liner  collapse  and  projectile  formation 
Design  (a)  resulted  in  excessive  stretching  of  the  peripheral  penetrator,  which  led  to  formation  of  a 
conoidal  debris  cloud  instead  of  the  desired  robust  penetrator.  Design  (b),  which  produced  a  robust 
hole  cleaner  and  stretching  hole  driller,  is  quite  satisfactory 
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Figure  5 

Liner  peripheral  portion  can  be  represented  as  an  idealized  two-mass  system  with  masses  m]  and 

m2  placed  at  the  edges  1  and  2 


Figure  6 

Liner  profile  can  be  changed  by  small  changes  in  edge  axial  position  Szx  and  &2  (a)  arid  in  a  liner 
edge  angle  a\  (b),  while  the  radial  positions  r\  and  r2  of  the  edges  are  fixed 
Small  adjustments  Sa \  in  liner  edge  angle  a\  result  in  negligible  radial  changes 

Sr  of  the  liner  profile. 
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Design 

CXl 

02 

<i 

h 

degrees 

cm 

A 

80 

74 

0.388 

0.388 

B 

80 

69 

0.388 

0.388 

C 

81 

63 

0.388 

0.349 

D 

81 

63 

0.370 

0.330 

E 

87 

58 

0.388 

0.349 

F 

81 

63 

0.370 

0.330 

G 

87 

58 

0.218 

0.333 

Calculated  hole  cleaner  shape  changes  due  to  changes  in  liner  (peripheral  portion)  profile 
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Figure  8 

Schematic  of  the  experimental  setup 
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Figure  9 

Radiographic  images  of  the  projectile  show  that  the  iiner  formed  a  robust  hole  cleaner 

and  a  stretching  hole  driller 


Hole  vroduced  by  ” Cleaner ” 


Hole  vroduced  by  Dnlier 


Figure  10 

Photograph  of  a  schematic  of  the  hole  produced  in  a  stack  of  steel  witness  plates 
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